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ABSTRACT 

Ultrasonic  waves  provide  an  excellent  means  to 
study  the  mechanical  properties  of  both  heterogeneous 
and  composite  materials.  Recently,  both  theoretically 
and  experimentally,  we  studied  characteristics  of 
elastic  wave  propagation  through  fiber  and  particle 
reinforced  composites.  Theoretically  we  analyzed  phase 
velocities  of  longitudinal  and  shear  waves  propagating 
through  composite  media  with  aligned  continuous  fibers 
and  ellipsoidal  shaped  particles.  Alignment  of  rein* 
forcing  fibers  (or  particles)  Imparts  anisotropy  to 
the  effective  mechanical  and  physical  properties  of  the 
composite.  Theoretically,  we  modeled  the  anisotropic 
phase  velocities,  which  agree  with  experimental 
measurements.  One  important  feature  of  this  combined 
theoretical  and  experimental  study  is  the  possibility 
of  obtaining  the  (often  uncertain)  fiber  (or  particle) 
elastic  properties  by  comparing  the  model  predictions 
with  the  observations. 

Ue  show  that  one  can  use  the  theoretical  model  to 
obtain  effective  properties  of  media  containing  voids 
(pores).  As  examples,  we  consider  porous  olivine  and 
creep  cavities  in  copper.  In  the  latter  case  some  ex¬ 
perimental  results  are  presented  which  show  good  agree¬ 
ment  with  predictions. 

Finally,  we  present  some  results  for  attenuation 
in  a  particulate  composite  Including  the  effect  of  non¬ 
ideal  Interface  properties.  . 

INTRODUCTION  ^  ~~ 

Determination  of  effective  elastic  moduli  and 
damping  properties  of  a  heterogeneous  material  by  using 
elastic  waves  (propagating  or  standing)  is  very  effec¬ 
tive.  Several  theoretical  studies  show  that  for  long 
wavelengths  one  can  calculate  the  effective  wave  speeds 
of  plane  longitudinal  and  shear  waves  through  a  com¬ 
posite  material.  At  long  wavelengths  the  wave  speeds 
thus  calculated  are  non-dlspersive  and  hence  provide 
the  values  for  the  static  effective  elastic  properties. 
References  to  some  of  the  recent  theoretical  and  exper¬ 
imental  works  can  be  found  in  (1-13).  The  scattering 
formulations  developed  in  (l-8)~provide  a  means  to 
obtain  not  only  the  effect! veTwave  speeds  but  also  the 
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damping  of  wave  amplitudes  due  to  scattering. 

In  this  paper  we  present  results  of  some  of  our 
recent  investigations  of  phase  velocity  and  attenuation 
of  plane  longitudinal  and  shear  waves  propagating  in  a 
medium  with  microstructure.  Microstructures  studied 
were  either  inclusions  or  fibers.  In  the  former  case, 
we  examined  the  effect  of  Inclusion  shape,  orientation, 
volume  fraction,  and  elastic  properties  on  wave  speeds. 
For  fiber-reinforced  materials  we  studied  continuous 
aligned  fibers.  In  either  case  the  medium  behaves 
anisotropically  because  of  the  alignment  of  the  inclu¬ 
sions  or  the  fibers. 

The  theoretical  model  used  a  wave-scattering 
approach  together  with  Lax’s  quasi-crystalline  approxi¬ 
mation  and  predicted  the  macroscopic  isotropic  elastic 
properties  for  the  case  of  random  orientation  of  Inclu¬ 
sions  and  anisotropic  elastic  properties  caused  by 
preferred  orientation.  Both  homogeneous  and  non-horo- 
geneous  distributions  of  inclusions  or  fibers  were 
considered.  The  scattering  approach  led  also  to  an 
estimation  of  attenuation  via  the  use  of  optical 
theorem. 

TFe  experimental  methods  consisted  of  a  pulse-echo 
technique  and  the  resonance  method.  These  were  chosen 
to  provide  the  advantages  of  small  specimens  and  low 
Inaccuracy.  For  details  of  the  experimental  techniques 
the  reader  is  referred  to  ( 10-.1£,15_). 

Both  homogeneous  and  nonhomogeneous  distributions 
of  Inclusions  and  fibers  were  considered.  First,  we 
examined  the  case  of  a  random  homogeneous  distribution 
of  randomly  oriented  spheroidal  inclusions  in  a  homo¬ 
geneous  matrix.  Second,  we  considered  a  random  homo¬ 
geneous  distribution  of  oriented  spheroids.  In  the 
first  case  the  macroscopic  properties  of  the  composite 
are  Isotropic.  In  the  second  aligned  spheroidal  inclu¬ 
sions  Impart  anisotropy  to  the  composite.  In  order  to 
model  the  SIC-particle-reinforced  composite  shown  in 
Fig.  1  we  used  a  two-step  process  combining  the  two 
distributions  described  above.  In  this  example  the 
material  consists  of  "Islands”  of  oriented  oblate 
spheroidal  A1  inclusions  in  a  "sea"  of  randomly  oriented 
SIC  prolate  spheroidal  Inclusions.  We  also  examined  by 
the  two-step  process  mentioned  above  the  fiber-rein¬ 
forced  composite  shown  In  Fig.  2  which  consists  of  a 
nonhomogeneous  distribution  of  aligned  continuous 


1  0  0  Jim 


A1 2®3  M  matrix. 

In  the  following  we  describe  the  theoretical  tech* 
nique  first  in  connection  with  a  composite  consisting 
of  ellipsoidal  shaped  inclusions  in  a  homogeneous 
matrix.  Special ization  to  the  two-dimensional  case  of 
fiber-reinforced  materials  is  then  briefly  discussed. 
Model  predictions  are  then  compared  with  experimental 
results,  finally,  we  present  some  model  calculations 
for  attenuation  in  a  particle-reinforced  composite  a-d 
for  elastic  properties  of  materials  with  voids. 

MULTIPLE  SCATTERING  BY  A  DISTRIBUTION  OP  ELLIPSOIDAL 
INCLUSIONS 

In  this  section  the  scattered  field  at  any  point 
in  the  matrix  is  obtained  in  the  presence  of  a  distri¬ 
bution  of  N  ellipsoidal  inclusions.  We  consider  bof. 
aligned  and  nonaligned  cases.  The  results  obtained 
are  valid  when  the  wavelength  is  long  compared  to  the 
Inclusion  dimensions.  The  method  represents  an  exte-- 
sion  to  the  elastic  case  of  the  approach  used  in  (16. 
17).  - 

In  expressing  the  scattered  field  in  the  present* 
of  a  number  of  ellipsoidal  inclusions  we  must  recall 
the  results  for  a  single  ellipsoidal  inclusion  (18). 

Let  the  center  of  the  i-th  ellipsoid  be  located  at 
(Xf,  Yj,  Zi)  referred  to  a  Cartesian  frame  of  refere'te 
axes;  let  the  principal  axes  of  this  ellipsoid  be 
obtained  from  the  Cartesian  axes  (XYZ)  by  rotation 
defined  by  the  Eulerian  angles  a; ,  ftj ,  yj.  Let 
u^(Plei)  denote  the  total  field  incident  on  the  i-th 
inclusion.  In  terms  of  vector  spherical  wave  functions 
that  are  regular  at  ej,  this  “exciting"  field  may  be 
written  quite  generally  as 


Fig.  1  Photomicrograph  of  SiC/At-al loy  composite 

formed  from  powders  of  SIC  and  At  alloy.  For 
present  purposes  we  represent  the  SIC  particles 
as  prolate  ellipsoids  with  an  aspect  ratio  of 
3.0.  Ue  represent  the  large  At  "islands”  as 
oblate  spheroids  (aligned  in  the  rolling  plane) 
with  an  aspect  ratio  of  0.33.  The  randomly 
oriented  SIC  particles  exist  in  a  "sea"  where 
the  SIC  volume  fraction  Is  approximately  double 
that  of  the  macroscopic  average  volume  fraction. 
The  SIC  particles  range  up  to  5  vm  in  their 
largest  dimension.  The  At  "island"  particles 
range  up  to  60  urn. 


Photomicrograph  of  Al,03Mt  alloy  composite. 
Uniaxial  Al^Oj  fibers  are  20  •  S  wm  in  dia- 

Note  nonhomogeneous  fiber  distribution. 
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where  rj,  ei,  are  the  spherical  polar  coordinates 
of  the  point  P  with 

V  t(X  -  C,)Z  ♦  (Y  -  ni)z  ♦  (Z  -  Cl)Z) 

Z  -  C1  «  rfCOSOj,  X  -  •  r<sinOicoso1 , 

Y  -  nj  *  r1sinojsin«i  (E) 

position  vector  R  referred  to  the  center  of  the 
ellipsoid.  Thus,  in  writing  Eq.  (l)  we  assume  that  the 
time  dependence  is  through  the  factor  e-’-'t,  which  hes 
been  dropped.  The  vector  wave  functions  appearing  ir. 
Eq.  (1)  were  given  in  (18).  t  is  the  ratio  of  the 
longitudinal -wave  and  sfiear-wave  speeds  in  the  matrix. 
The  field  given  by  Eq.  (1)  will  be  scattered  by  the 
ellipsoid  and  this  scattered  field,  denoted  by 


"h«re  ^iuv  and  -iuv  are  th*  appropriate  spherical 

vector  wave  functions  that  satisfy  the  radiation  condi¬ 
tions  as  rt—.  Here  c  •  wc/cj,  where  cj  denotes  the 


u(3) 


meter. 


longitudinal  wave  speed  in  the'matrix  and  the  principal 
semi  axes  of  the  ellipsoid  are,a,b,c(a>b»c).  Equations 
(1)  and  (3)  keep  terms  to  0{e3). 

The  relationships  between  the  scattered*  field 
coefficients  A^v,  8<uv  and  the  incident*fie1d  coeffi¬ 
cients  ajm,  b^,  and  C{m  were  derived  in  (jg)  when 
the  semiaxes  of  the  ellipsoid  were  parallel  to  the 
Cartesian  axes  XYZ.  For  aligned  ellipsoids  this  choice 
can  be  made  without  loss  of  generality.  However,  for  a 
nonaligned  configuration  it  is  necessary  to  obtain 
these  relationships  for  arbitrary  orientation.  To 
derive  these  general  relations  from  the  results  of  Qfi) 
it  is  necessary  to  refer  both  the  incident  field  (1) 
and  the  scattered  field  (2)  to  the  ellipsoid  axes. 

This  is  done  by  using  the  rotational  addition  theorems 
for  spherical  harmonics  (19.201.  One  can  show  that 

uE(R|Rf)  *  I  l  [«.  ill!.  (r.,o;.o;)  ♦ 

Tbim‘n  ?im’n  <ri  ,ei  }  * 

TCim'n  -im'n  (ri,0i'*i}  J  <4) 

In  writing  Eq.  (4)  the  relations  expressing  the  spheri¬ 
cal  vector  wave  functions  referred  to  XYZ  and  the 
ellipsoidal  axes  have  been  used.  These  have  the  form: 

(1)  n 

Liwi  *  l  B(m,m*,n)  L  .Jr.e'.e')  (5) 

in' *-n 

Similar  expressions  hold  for  and  N^.  Here 

<fc3iK)i>1,z  «£?i  (*-)  w 

Using  Eq.  (5)  in  Eq.  (1)  it  follows  that 

41m'n  *  l  «imn  B(m.m-.n).  etc.  (7) 

in*«n 

The  scattered  field  referred  to  the  new  axes  is  then 


iv0c-  uri 

A,  •  '  ■  ■<  ?  I  T  .  (a,  ♦  «(v)b.  )  (ll) 

v  Aiic3  vj  H  w  ‘Vl  Vl 


4(v)  *  3r  ,  v  »  0,2 
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Explicit  expressions  for  T“l^l  are  given  in  (9). 

Using  Eqs.  (7)  and  (11)  in  Eq.  (10)  it  Is  found 

that 
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(12) 
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This  completes  the  derivation  of  the  scattered  field 
due  to  a  single  arbitrarily  oriented  ellipsoid. 

He  now  seek  the  scattered  field  due  to  a  number  of 
arbitrarily  oriented  ellipsoids.  It  is  necessary  to 
expand  the  scattered  field  due  to  the  j-th  ellipsoid 
in  vector  wave  functions  regular  near  the  i-th  one. 
Using  translational  addition  theorems  (21)  the  exciting 
field  on  the  1-th  ellipsoid  is 

uE(R|Ci)  ■  u(1)  ♦  f  J  [a2  l(l)  ♦ 
n*0  m*-n  imn  *imn 

tb?  1)  ♦  TC^  m(D,  1 1 

Dimn  .imn  cimn  .imn* 
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and  the  asterisk  denotes  complex  conjugate.  Thus, 
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Here  u(1)  is  the  Incident  field  and  Au‘w  and  Buv  are 
•  mn  mn 

defined  in  (2). 

The  incident  wave  will  be  taken  as  a  plane  wave 
propagating  along  the  Z*axis.  Thus,  it  may  be  written 

as 

(1)  1klZ  1k2Z 
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Using  Eqs.  (IS)  and  (16)  in  Eq.  (12)  the  scattered- 
field  coefficients  correct  to  0(c3)  due  to  a  number  of 
arbitrarily  oriented  ellipsoids. are  obtained.  In  this 
study  the  calculations  are  simplified  If  we  assume  that 
the  distribution  of  particles  is  random  and  homogeneous. 
Me  assume  further  that  the  particles  have  the  same 
shape  and  size. 

Frequently,  the  processing  of  these  composites 
produces  a  material  that  contains  non  homogeneities,  for 
example  matrix  Islands  devoid  of  particles  (see  Fig.  1). 
This  fact  Is  considered  In  a  second-step  calculation 
that  assumes  that  there  Is  a  distribution  of  aligned 
oblate-spheroidal  matrix-material  Inclusions  In  the 
composite. 

EFFECTIVE  WAVE  SPEEDS  IN  A  COMPOSITE  OF  RANDOM 
DISTRIBUTION  OF  DISORIENTED  ELLIPSOIDS 

To  calculate  the  effective  wave  speeds  of  plane 
waves,  we  take  an  ensemble  average  of  Eq.  (12).  For 
this  purpose,  the  probability  density  of  finding  the 
scatterers  at  Rj,  Rj,  ...  Is  denoted  by  p(R^,...,RN), 

which  can  be  written 
p(Rj»...»  Rn)  *  P  (H|)  p(R2>. •  • » 

•  P(R1)  P<R2IR1)  p  <«3 . V?!-  ?2}  (173 

Assuming  the  distribution  to  be  uniform  It  follows  that 

p(Rj)  •  1/V  ,  RjtV 

•  0  ,  Rj<V  (18) 

where  V  is  the  volume  of  the  composite.  Furthermore, 
following  the  analyses  of  (3),  the  conditional  density 
p(Rj|Rj)  will  be  taken 

P(«2J«1)  •  1/V  ,  |R2  -  Rj|  >  2a 

•  °  .  |Rj  *  Rjl  <  2a  (19) 

This  it  the  "well -stirred"  approximation  whose  validity 
for  high  concentrations  Is  open  to  question  (7).  How¬ 
ever,  the  results  obtained  In  (£)  under  this  assumption 


coincide  with  the  lower  (upper)  bound  of  Hashln  and 
Shtrlkman  (22).  (See  Ref.  (3).)  Although  this  is  no 
justification  for  adopting  Eq.  (19),  we  use  It  here  for 
slmpl iclty. 

In  taking  the  ensemble  average  cf  Eq.  (12),  for 

U«V« 

random  orientation,  the  average  of  Tu*  1  over  all 
orientations  Is  simply 

ulv.  v  , 

«T  1  *>  •  3  p  •  l  T*TT  T  1  (20) 
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Here,  ng  Is  the  number  density  of  the  SiC  particles. 

Equation  (21)  has  the  same  form  as  derived  in  (2) 
for  a  distribution  of  spherical  particles.  It  contains 
the  conditional  expectations  with  two  particles  held 
fixed.  An  approximation  to  «Ajyv»j  is  obtained  by 

assuming  the  Lax  (£1)  quasi -crystal  1 ine  approximation. 
This  approximation  assumes  that 


<A,  *  «A .  >, 
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Substitution  of  Eq.  (22)  Into  Eq.  (21)  leads  to  an 
Integral  equation  for  the  determination  of  <Aj  »j.  To 

find  the  effective  wave  speeds  of  plane-wave  propaga¬ 
tion,  we  assume  a  solution  to  the  Integral  equation: 


<A1uv>1  " 
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uv 


where  K  is  the  effective  wavenumber  and  X  are 
constants. 

Substituting  Eq.  (23)  Into  Eq.  (21)  it  can  be 
shown  following  (2)  that  for  plane-longitudinal  and 
plane-shear  waves'the  respective  wavenumbers  K.  and  K. 
are  1  L 

K7  (  1+9cPj)  (l*3cP0)  (1  ♦  |  cP2  (2+3k|/k?)  1 

~7  *  - 1 - : — - r~z - 5 — 5 —  (24) 

kj  1-15cP2  (1*3cP0)  ♦  j  cP2  (2*3k‘/k‘) 


(  1+9cPj)  (  1  +  |  cP2  (2+3k|/k|) ) 
kj  1  ♦  |  CP2  (4-9k*/k*) 


c  ■  y  «abcnQ. 

ull 

Using  the  expressions  for  T  1  (9)  It  Is  found  that 

“l1 

9Pj  •  (o’/o-  1)  (26) 

Now  writing 

Kj  *  t/Cj  •  K2  *  w/C£ 

where  Cj*.  c2*  are  the  effective  longitudinal -wave  and 

and  shear-wave  speeds,  the  effective  Lame  constants  x* 
and  u*  are  obtained  from  the  equations 

>1.  .  <1**V  11  *  I  lf2 

x  Zu  1-15cP2  (1*3cP0)  *  \  cP2(2*3kz/kz) 

*  1  ♦  l  cP,  (2+3k?/k?) 

•  — hr2 - r4*  (28) 

“  1  *  f  cP2  (4-9k|/k‘) 

The  expressions  given  a^ove  can  be  simplified  by 
using  the  expressions  for  Tl'lW  In  (20).  The  effective 
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bulk  modulus  and  shear  modulus  are 
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Taking  now  the  ensemble  average  of  Eq.  (12)  one  finds 
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As  before,  the  quasi-crystalline  approximation  is 
made  to  derive  an  integral  equation  for  Consi¬ 

dering  now  a  longitudinal -plane-wave  solution  of  this 
equation  (<Ajgw»i  »  X^e1  Klci ;  u  •  0,  v  «  0,  1,  2)  a 
set  of  linear  homogeneous  equations  in  XQq,  XQj,  and 
XQ2  Is  obtained.  Keeping  only  the  lowest-order  terms 
2  2 

In  Kj/kj  and  equating  the  determinant  of  the  coeffi¬ 
cients  of  Xqq,  Xqj,  and  Xq2  to  zero  one  finds 


4“<l*3cT°°)  ((1*cT°°)  (1*cT°Z(1* 
*1 


u  -  li 
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The  tensor  T<jkl  has  been  defined  In  (9).  It  is  noted 

that  in  the  low- volume  concentration  limit  these  agree 
with  the  expressions  obtained  by  Boucher  (|4). 

EFFECTIVE  WAVE  SPEEOS  IN  A  COMPOSITE  OF  RANDOM 
DISTRIBUTION  OF  ORIENTED  EUIPSOIOS 

When  the  ellipsoids  are  aligned,  effective  wave 
speeds  of  plane  waves  can  be  calculated  in  the  same 
manner  as  In  the  previous  section  by  taking  the 
enseeble  average  of  Eq.  (12).  Hdwever,  for  propagation 
in  an  arbitrary  direction  the  dispersion  equation  Is 
complicated.  The  equation  simplifies  for  propagation 
along  the  ellipsoidal  axes.  Results  for  this  particu¬ 
lar  case  are  given  here.  Details  of  the  derivation  are 
omitted. 

When  an  incident  plane  wave  propagates  along  one 
of  the  ellipsoid  axes,  c-axis  say,  the  XT2  axes  can  be 
chosen  tobe  parallel  to  the  elllosoidal  abc-axes.  Then 

•ui v.  ujv, 

T  1  1  •  T  1  1  Mil 
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♦  T02  +  5  S  TSo  TS?>J  (33) 

For  a  plane  shear  wave  polarized  in  the  X-direc* 
tion  and  propagating  along  the  Z-direction  the  disper¬ 
sion  equation  is  obtained  by  taking 

1K»c( 
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It  Is  found  that 
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It  is  easily  shown  that  for  a  shear  wave  polarized 
along  the  Y-axis  and  propagating  along  the  Z-axis  the 
effective  wavenunber  k2  Is  obtained  by  replacing 

T12  *  6  T*122  w1th  T12  '  6T*12?  1n  Eq-  (34)< 

The  above  derivation  Is  for  propagation  along  the 
Z-axis.  Results  for  propagation  along  the  X-axis  and 


Y-axis  are  obtained  by  cyclic  interchange  of  abc.  It 
Is  shown  easily  that  the  speed  of  propagation  of  a 
shear  wave  polarized  along  the  X-direction  and  moving 
along  the  2-direction  is  the  same  as  that  of  a  shear 
wave  polarized  along  the  2-direction  and  moving  along 
the  X-direction. 

WAVE  PROPAGATION  IN  A  FIBER-REINFORCED  COMPOSITE 

The  analysis  given  above  can  be  modified  easily  to 
treat  wave  propagation  in  a  plane  perpendicular  to  the 
fibers.  Propagation  of  longitudinal  and  shear  waves 
polarized  and  propagating  in  this  plane  was  considered 
in  (1).  Although  the  treatment  there  was  for  isotropic 
fibers,  the  extension  to  transversely  isotropic  fibers 
is  made  easily.  This  was  done  in  (2i).  In  that  study 
the  problem  of  a  shear  wave  polarized  along  the  fibers 
and  propagating  perpendicular  to  them  was  also  con¬ 
sidered. 

Talcing  the  2-axis  along  the  fiber  axes  and 
denoting  the  five  independent  elastic  constants  charac¬ 
terizing  the  fibers  as  C^,  Cj2>  C13«  C33>  and  C44,  it 

was  shown  (2£)  for  SH-wave  propagation  that  the  effec¬ 
tive  wave  number  8*  is 


EXPERIMENT 

Si  11 con-Carbide/ Aluminum 

The  photomicrograph  shown  in  Fig .  1  represents  an 
SiC-  re  in  forced  A?  composite.  The  material  was  obtained 
from  a  commercial  supplier  in  the  form  of  1-cm  plates. 
Nominally,  the  plate  contained  30  volume  percent  SiC. 

Sound  velocities  were  determined  by  a  pulse-echo 
method  described  in  detail  previously  (li).  Briefly, 
1-cm  cubes  were  prepared  by  grinding  so  that  opposite 
faces  were  flat  and  parallel  within  5  urn.  Quartz 
piezoelectric  crystals  with  fundamental  resonances 
between  4  and  7  MHz  were  cemented  with  phenyl  sali¬ 
cylate  to  the  specimens.  An  x-cut  transducer  was  used 
for  longitudinal  waves  and  an  ac-cut  for  transverse 
waves.  Ultrasonic  pulses  1  to  2  cycles  long  were 
launched  into  the  specimen  by  electrically  exciting  the 
transducer.  The  pulses  propagated  through  the  speci¬ 
men,  reflected  from  the  opposite  face,  and  propagated 
back  and  forth.  The  pulse  echoes  were  detected  by  the 
transducer  and  displayed  on  an  oscilloscope  equipped 
with  a  time  delay  and  a  microprocessor  for  time- 
interval  measurements.  The  sound  velocity  was  computed 
by 


where  m  «  C44/u,  6*  *  j) ,  o*  the  effective 

density,  and  uLT  the  effective  axial  shear  modulus. 

For  propagation  of  longitudinal  and  shear  waves  polar¬ 
ized  and  propagating  in  a  plane  perpendicular  to  the 
fibers  it  was  found  that  the  corresponding  wavenumbers 
are 

K1  aMl  ♦  c  P  )J1  ♦  e  Pja+kf/kJ)] 
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where 

kt  -  (*+u)  u(Cfi»  -u) 

P  S  .  — 1 _  '  p  ■  .  66 _ 

S  *T  *  u  Cgg(x+3u)  ♦  u(x+u) 

Kj  is  the  plane-strain  bulk  modulus  of  the  fiber 

Although  the  above  derivation  is  for  transversely 
isotropic  fibers  in  an  isotropic  matrix.  It  is  easily 
generalized  to  the  case  of  Isotropic  fibers  in  a  trans¬ 
versely  isotropic  matrix  with  the  symmetry  axes 
parallel  to  the  fibers.  This  1$  needed  to  deal  with 
the  situation  shown  in  Fig.  2.  As  in  the  case  of  the 
inclusion  problem  described  above,  it  is  necessary 
first  to  calculate  the  moduli  of  the  effective  trans¬ 
versely  Isotropic  composite  medium  formed  by  a  random 
homogeneous  distribution  of  fibers  in  a  homogeneous 
isotropic  matrix.  In  the  second  step  we  Introduced  a 
distribution  of  homogeneous  isotropic  matrix  "fibers" 
into  the  composite  constructed  in  the  first  step. 


v  «  2t/t 

where  i  denotes  specimen  length,  and  t  the  round-trip 
transit  time.  On  the  oscilloscope,  t  was  the  time 
between  adjacent  echoes,  the  first  and  second  echoes 
usually  being  measured,  and  within  these  the  time 
between  leading  cycles.  Elastic  constants  were 
computed  from  the  general  relationship 

C  »  ov2 

where  p  denotes  mass  density. 

Mass  density  was  determined  by  Archimedes's  method 
using  distilled  water  as  a  standard.  For  30  volume- 
percent  SiC/Al,  we  found  o  *  2.838  g/enn. 

Using  quantitative-metallographic  equipment  we 
verified  the  overall  SIC  volume  fraction  to  be 

c  •  28.1  i  1.9  percent  and  the  Ai-island  concentration 
to  be  c*  m  48.5  percent.  This  means  that  within  the 
sea  the  SIC  volume  fraction  was  54.6  percent.  For  the 

calculations  reported  below,  we  took  c  *  0.30  and 
varied  c*  from  0.0  to  0.5.  Metallography  showed  that 
the  SiC  particles  represented  as  prolate  spheroids 
possess  an  asepet  ratio  of  approximately  3.0.  Simi¬ 
larly,  the  aluminum  islands  represented  as  oblate 
spheroids  possess  an  aspect  ratio  of  approximately 
0.33. 

Sapphire  (A».203)/Alum1num  and  Boron/Aluminum 

Figure  2  shows  the  microstructure  of  sapphi re 
fiber-reinforced  aluminum  material,  which  was  obtained 
from  a  cofimerlcal  supplier  as  6-nm  plates.  The  matrix 
consists  of  a  2XLi-At  alloy.  The  fibers,  553  by 
volume,  consist  of  99+3  o-alumina,  20+5  um  diameter, 
with  a  manufacturer-reported  mass  density  of  3.95  g/cm3. 
Sound  velocities,  volume  fraction,  and  mass  density  of 
the  composite  were  determined  as  described  above  for 
SIC/At. 

Elastic  constants  of  the  fiber  were  estimated  from 
Tefft’s  (2£)  monocrystal  trigonal -symmetry  elastic 
constants.  These  were  averaged  to  the  quasi-isotropic 
polycrystal  case  by  a  Voigt-Reuss-Hill  arithmetic 
average.  The  polycrystalline  values  were  then  scaled 
downward  slightly  to  agree  with  the  observed  S33,  the 

reciprocal  Young's  modulus  along  the  fiber  direction. 

For  this  scaling  down,  Poisson's  ratio  was  kept  fixed. 


This  gave  a  fiber  Young's  modulus  of  3S8  GPa,  within  the 
340*380  G?a  range  specified  by  the  fiber  manufacturer. 
This  estimate  assumes  isotropic  elastic  properties  in 
the  fiber.  In  addition  to  the  pul$e*echo  method  we 
also  used  the  resonance  method  with  a  three-component 
(Marx)  oscillator  to  measure  the  elastic  properties  of 
Boron-Aluminum  composite.  The  details  of  this  can  be 
found  in  (27). 

Creep-Cavities  in  Copper 

Ultrasonic  pulse-echo  method  was  used  to  measure 
longitudinal  sound-wave  velocities  along  and  perpendic¬ 
ular  to  the  axis  of  stress  applied  to  produce  creep  in 
polycrystalline  copper.  The  creep  cavitation  produced 
is  shown  in  Fig.  3. 

RESULTS:  CALCULATIONS  ANO  OBSERVATIONS 

Ef fecti ve  El asti c  Properti es 

For  SiC/Ai,  table  1  shows  the  principal  results 
of  the  study  represented  as  elastic-stiffness  coeffi¬ 
cients,  Cjj  in  Voigt's  notation.  Columns  2  and  3  give 

the  properties  of  the  constituents:  aluminum  alloy 
6061  (measured  in  this  study)  and  SIC  (reported  by 
Schreiber  and  Soga  (26)  ).  Column  4  contains  the 
observed  elastic  constants  of  30- volume-percent  SiC/AI. 
Column  6  gives  the  prediction  based  on  the  assumption 
of  homogeneously  distributed  SIC  particles.  Finally, 
column  7  gives  predictions  based  on  the  present  model 
for  the  nonhomogeneous  case  shown  in  Fig.  1.  The 
co-ordinate  system  is  x3  perpendicular  to  plate,  In 

plate  in  rolling  direction,  x?  in  plate  perpendicular 
to  rolling  direction. 


Table  2  -  Observed  and  Predicted  Static  Elastic 

Stiffness,  C^j,  of  55-volume-percent  Sapphire 

Fiber/Aluminum  at  Ambient  Temperature. 


A* 

— — — r 
Sapphire1 

i 

Obs. 

Homog. 

Nonhomog. 

Spec. 

grav. 

2.700 

3.95 

3.22 

3.39 

3.39 

C!1 

1.094 

4.193  | 

1.92 

1.971 

•  1.985 

C12 

0.560 

1.273  I 

0.66 

0.855 

0.857 

C13 

0.560 

1.273  1 

0.67 

0.780 

0.782 

C33 

1.094 

4.193  , 

2.65 

2.730 

2.730 

^44 

0.267 

1.46 

0.64 

0.593 

0.603 

C66 

0.267 

1.46  | 

0.61 

0.558 

_ 

0.564 

Predictions  in  this  table  are  based  on  c  *  0.55  and 
c*  *  0.10.  Units  are  lO^N/m2. 


Table  2  gives  similar  results  for  the  Al^Oj/A: 

fiber-reinforced  composite.  The  co-ordinate  system  is 
Xj  in  fiber  direction.  For  comparison  with  theory,  we 

“forced"  a  transverse-isotropic  symmetry  on  the  obser¬ 
vations  by  the  relationships  Cj 3  «  (Cjj  ♦  C22)obs  ^ ' 

C!3  *  <C13  *  C?3W'2-  C«4  '  <C«  •  CSsW/2 
C66  "  (2C66  4  C11  ‘  C12lobs./*' 


Table  I  -  Observed  and  Predicted  Elastic  Stiffnesses,  C^,  of  30- Volume-Percent  Si  C/606 1  A1  at  Room  Temperature. 


6061 

SIC4 

30- vol  ume -percent  SiC/6061  A1 

Obs. 

R-O-fll 

Homog. 

Nonhomog. 

Spec. 

Grav. 

2.706 

3.181 

2.838 

2.849 

2.849 

2.849 

C11 

1.105 

4.742 

1.659 

2.195 

1.583 

1.727 

C22 

1.105 

4.742 

1.651 

2.195 

1.583 

1.727 

C33 

1.105 

4.742 

1.483 

2.195 

1.583 

1.480 

C44 

0.267 

1.881 

0.433 

0.751 

0.443 

0.429 

C55 

0.267 

1.881 

0.487 

0.751 

0.443 

0.429 

C66 

0.267 

1.881 

0.487 

0.751 

0.443 

0.513 

C12 

0.571 

0.980 

0.685 

0.693 

0.697 

0.701 

C13 

0.571 

0.980 

0.662 

0.693 

0.697 

0.677 

r23 

0.571 

0.980 

0.662 

0.693 

0.697 

0.677c 

4  From  Schreiber  and  Soga  (28) . 
6  He  assume  C23  *  Cj3. 


c  We  assume  C13  •  C23  * 


12’ 


Predictions  in  this  table  are  based  on  c  «  0.30,  c*  •  0.50, 
SiC-particle  aspect  ratio  *  3.0,  At-island  aspect  ratio  •  0.33. 

Units  on  are  10llN/m2. 


0.966  C 


for  a  fiber-reinforced  8oron-Aluminum  results  of 
observation  and  predictions  of  the  random-distribution 
model  used  in  (J)  and  two  other  periodic  distribution 
models  are  shown  in  Table  3.  In  the  model  calculations 
both  constituents  were  assumed  to  be  isotropic,  for 
the  Aluminum  matrix  this  assumption  was  verified  exper¬ 
imentally.  However,  for  Boron  fibers  the  elastic  con¬ 
stants  are  less  certain,  a  wide  range  of  values  being 
reported  in  the  literature.  The  values  shown  in  Table 
3  arose  from  fitting  a  linear  rule-of-mixtures  to  our 
observed  value  combined  with  Gschneider's  (22) 

recommended  values  of  Boron's  bulk  modulus.  It  may  be 
noted  that  one  can  derive  the  properties  of  Boron  by 
using  the  observed  values  and  the  predictions  of  the 
random  distribution  model.  If  that  is  done  then  It  is 
found  that  Boron  fibers  are  slightly  anisotropic  with 

11  2 

the  properties  (all  in  units  of  10  N/m  ): 
kT  »  2.693.  uLT  -  1.891,  un  *  2.014 

El  *  3.72,  vLT  *  0.13,  vT  -  0.121,  Et  •  4.515 


Young's  modulus  of  cast  iron  has  been  studied  experi- 
mentally 'in  (32,).  Similar  studies  were  also  reported 
in  (34.351.  Fig.  5  shows  this  dependence  using  our 
model.  Also  shown  in  this  figure  are  the  results  of 
various  experimental  studies.  We  made  two  model  calcu¬ 
lations  corresponding  to  the  properties  of  graphite  by 
the  lower  and  upper  bounds  (22).  As  seen  from  this 
figure  lower  bound  properties  of  graphite  give  results 
closer  to  the  experimental  observations. 

Attenuation 

In  this  section  we  derive  expressions  for  the 
attenuation  coefficients  in  a  medium  containing  a 
distribution  of  spherical  inclusions  with  thin  inter¬ 
face  layers  through  which  the  elastic  properties  vary 
rapidly  from  those  of  of  the  inclusions  to  those  of  the 
matrix.  Such  interface  layers  are  often  present  due  to 
processing  (see  for  example  (37,32)  )• 

Consider  a  distribution  of  spherical  inclusions  of 
elastic  properties  Xj,  uj,  and  density  Cj  enbedded  in  a 

matrix  material  with  properties  x^,  and  o,.  Let 


Table  3  •  Observed  and  Predicted  Static  Elastic  Properties  of  48X  Boron  Fiber/ Aluminum  Composite  at  Ambient 
Temperature 


At 

Boron 

Observed 

Square 

Model4 

Hexagonal 
Model b 

Random 
Model c 

Full  Random 
Model 

C11 

1.107 

4.116 

1.852 

1.856 

1.872 

1.790 

1.790 

C12 

0.573 

0.590 

0.779 

0.661 

0.745 

C13 

0.573 

0.590 

0.606 

0.578 

0.583 

C33 

1.107 

4.116 

2.450 

2.480 

2.551 

2.560 

C44 

0.267 

1.763 

0.566 

0.451 

0.561 

0.559 

0.559 

C66 

0.267 

1.763 

0.526 

0.606 

0.523 

0.523 

*  After  Achenbach  (JO) 
b  After  Hlavacek  (Jl) 
c  After  (J) 

Wave  velocities  along  the  c-axis  of  oriented  ellipsoidal 
inclusions  as  given  by  Eqs.  (33)  and  (34)  are  dependent 
on  the  shapes  and  elastic  properties  of  the  ellipsoids. 
From  these  one  can  derive  the  results  for  a  medium  with 
oriented  ellipsoidal  shaped  cavities.  Such  is  the  case 
depicted  in  Fig.  3.  We  show  calculated  longitudinal 
wave  velocity  in  the  x^  direction  in  Fig.  4.  It  is 

assumed  that  the  void  aspect  ratio  is  1:10  (c/a). 

Results  for  spherical  voids  are  also  shown  together 
with  the  experimentally  measured  velocities  in  the  x^ 

(stress  axis)  and  Xj  directions.  This  figure  shows 

that  the  material  had  some  initial  anisotropy  due  to 
texture  which  has  not  been  taken  into  account  in  the 
modeling.  Assuming  completely  oriented  disc  shaped 
voids  one  can  estimate  their  aspect  ratio  by  determin¬ 
ing  which  matches  the  observed  curves.  For  v3  we 

obtain  1:7;  for  Vj  we  get  1:12;  thus  an  effective 

estimate  is  1:9.  More  details  of  the  experiment  can 
be  found  in  (J2). 

The  aspect  ratio  dependence  of  elastic  velocities 
in  particulate  coaposites  as  predicted  by  equations 
(24)  and  (25)  was  studied  further  in  the  context  of 
cast  iron.  Effect  of  graphite  particle  shape  on 


each  Inclusion  be  separated  from  the  mat,  x  by  a 
uniform  thin  layer  of  thickness  h  (<<a)  and  variable 
material  properties  x(r),  u(r),  and  o(r).  Here  r  is 
the  distance  from  the  center  of  the  sphere  and  a  is  its 
radius.  Because  of  the  synmetry  the  scattered  field  due 
to  a  single  sphere  can  be  calculated  in  an  exact  form. 
For  the  incident  field  given  by  (14)  it  can  be  shown 

that,  correct  to  0(c3),  the  scattered  field  is  given  by 
(3)  with 


X.  »  1cJ  (Pa.  ♦  Q  x.  ) 

iuv  1  V  1 UV  \  iuv 

< 

(no  sum  over  v) 


B,  „  *  1t2c3  (R  a.  ♦  S  X.  } 
iuv  v  iuv  v  iuv 


(38) 


where 


P  ■ 


1 


h  3X,  +  2u, 

3X,  ♦  2u-  -  (3X  +  2u. )  { 1  -  £  —6 - 6-  K.) 

i  c  1  1  a  +  2uJ  1 

dujMSXj^UjHl*  Kl) 


Fig.  3.  Photomicrograph  of  polycrystal  line  copper  con¬ 
taining  creep  cavities  at  grain  boundaries. 
Reference  mark  equals  0.5  nm.  Earlier  stage 
shows  small  disconnected  spherical  voids. 


Fig.  5.  Variation  of  Young  modulus  with  graphite- 

particle  aspect  ratio.  Symbols  represent  meas¬ 
urements.  Curves  represent  predictions  based 
on  present  model.  Predictions  occur  for  two 
volume  fractions:  ten  and  twelve  percent. 

Upper  curves  represent  upper  bounds  on  graph¬ 
ite's  elastic  stiffnesses;  lower  curves  repre¬ 
sent  lower  bounds. 


VOID  VOLUME  FRACTION,  c 


Fig.  4.  longitudinal  sound  velocity,  v,  versus  volume  fraction  of  voids,  c.  lines  1  and  3  represent 
velocities  in  directions  xj  and  *3,  where  X3  corresponds  to  the  tensile-stress  axis  and  x,  is 
perpendicular  to  X3.  Lines  A,  B,  and  C  represent  calculations:  A  for  spheres,  B  and  C  for 
oblate  spheroids  (discs)  with  a  1:10  aspect  ratio,  where  B  corresponds  to  randomly  oriented 
discs  and  C  to  discs  oriented  perpendicular  to  X3. 
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In  writing  the  expressions  for  and  we 

followed  the  notation  of  (1)  for  easy  comparison.  It 
was  assumed  In  (_1)  that  Kj~*  0  and  f  0.  In  the  case 

considered  here 

1  1 

K1  '  |  TTa^luT  ’  *2  *  |  g(a  ~  hx) 

where  we  assumed 

\{r)  ♦  2u(r)  ■  (\'x  ♦  2up  f( r)  ,  a  <.  r  <  a  +  h 
u(r)  •  u  *|  g(r)  ,  a  «  r  «  a  +  h 

Following  the  steps  leading  to  Eqs.  (24)  and  (25) 
It  is  then  found  that  the  effective  wave  speeds  are 
given  by  Eqs.  (24)  and  (25)  with  the  Py  defined  above. 

The  attenuation  caused  by  scattering  can  also  be 
calculated  by  noting  that  the  scattering  cross  sections 
for  P  and  S  waves  are 


p  4na4  1  « 


9v2  c4r 


(P  2  +  3(1  +  2t3)  P  2  ♦ 


5(1+|  tS)  P22] 


I  *  — —4- ■  (3Pj2  (1  +  2t3)  +  t2  P22  (1  ♦  |  r5)) 

4-a  (40) 

where  v„*  ^  a3.  The  attenuation  coefficients  a  and 
U  J  J 

a$  are  then  given  by 

“p  *  Vp  •  as  *  Vs-  (41) 

Note  that  attenuation  coefficients  derived  here 
are  for  lossless  inclusions  and  matrix.  If  there  is 
dissipation  in  either  of  these  materials  then  there 
will  be  additional  attenuation  which  may  be  derived 

from  (24)  and  (25)  by  assuming  that  \  and  u  are  complex 

(see,  for  example,  (40)  ). 

DISCUSSION  AND  CONCLUSION 

A  two-step  multiple-scattering  formalism  has  been 
used  to  predict  the  anisotropic  elastic  constants  of 
both  particle-reinforced  and  fiber- re  in  forced  composite 
materials.  Static  effective  properties  are  obtained  by 
taking  the  long-wavelength  limit.  For  both  materials, 
the  present  model,  involving  a  simple  nonhomogeneous 
distribution  of  particles  or  fibers,  both  explains  the 
elastic  anisotropy  and  gives  better  absolute  predic¬ 
tions  of  the  elastic  constants. 

In  addition  we  have  presented  results  for  phase 
velocities  in  a  fiber-reinforced  material  with 
anisotropic  fibers.  These  indicate  that  using  mooeiing 
and  observations  it  is  possible  to  infer  the  properties 
of  fibers  which  are  sometimes  hard  to  obtain.  This 
aspect  of  the  present  study  is  very  useful  in  other 
contexts  as  well,  like  cast  iron  with  graphite  parti¬ 
cles. 

The  modeling  reported  here  is  also  found  to  be 
applicable  to  medium  with  oriented  or  disoriented  voids. 

Finally,  we  have  presented  expressions  for  atten¬ 
uation  in  a  medium  reinforced  by  spherical  particles 
with  interface  layers. 
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